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Abstract. We study the asymmetric exclusion process on a regular Cayley tree with 
arbitrary co-ordination number. In this model particles can enter the system only 
at the parent site and exit from one of the sites at the last level. In the bulk they 
move from the occupied sites to one of their unoccupied downward neighbours, chosen 
randomly. We show that the steady state current that flow from one level to the 
next is independent of the exit rate, and increase monotonically with the entry rate 
and the co-ordination number. Unlike TASEP, the model has only one phase and 
the density profile show no boundary layers. We argue that in blood, air or water 
circulations systems branching is essential to maintain a free flow within the system 
which is independent of exit rates. 
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Exclusion processes pQ has been studied extensively as the paradigm models of non- 
equilibrium phase transitions. They exhibit reach variety of phases, phase-coexistence, 
shock-profiles and non-trivial boundary layers. Some variations of these models are 
exactly solvable [2] on a one dimensional lattice, which provide deep understanding of 
non-equilibrium transport, traffic and jamming. However, very little is known about the 
systems beyond one spatial dimension. In particular, transport in irregular structures, 
like networks has been a recent topic of interest [3J. In a generic undirected network 
particle can enter or exit at any arbitrary sites. Again, the presence of loops in the 
generic networks, also make the study of particle transport difficult. A prototype 
network is a Cayley tree, where the direction of transport, the entry and exit points are 
well defined. Absence of loops make the study relatively simpler. Again, several physical 
systems like, water transport in trees, transport of nutrients in blood-circulation system 
0], transport of antibody in idiotypic networks in immune system [5], air circulation in 
lung [6], and flow on disordered networks [7] are strikingly similar to this model system 
of Cayley trees. 




Figure 1. Rooted Cayley tree with q = 2. There are N levels and M — (2 N — 1) 
sites, labeled by i = 0, 1 . . . , N — 1 and j = 0, 1 . . . , M respectively. Particle enter at 
the parent site j = with rate a and exit from any of the sites at the last level N — 1 
with rate /3. 

In this Letter we study transport in a Cayley tree with arbitrary co-ordination 
number. Particles can enter the Cayley tree only at the parent site with rate a and 
are allowed to leave from any of the sites at the last level with rate (3. In the bulk, 
when allowed by hard-core interaction, the particle can move to one of its downward 
neighbour chosen randomly. The resulting current that flow between neighboring levels 
is found to be independent of the exit rate (3 when q > 2. Further, the current is found 
to be larger than that of the TASEP in one dimension for any value of (a, (3). The 
average density at each level, except the last one, do not show any finite size correction 
and decay exponentially as one moves away from the parent site. In the last level, the 
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density depends both on a and j3 and decays exponentially with the system size. We 
did a mean field analysis, which correctly reproduces all these features. 

First the model. Cayley tree, is a connected cycle-free graph. A iV-level Cayley 
tree, labeled as % = 0, . . . iV — 1, with co-ordination number z = q + 1 can be generated 
by starting form the parent level i = with one single site j = called root. The next 
generation of sites are then created iteratively, where each site in level i is connected to 
q new sites in the next level i + 1. Thus level i has q % sites and the resulting network 
has total M = (q N — l)/(q - 1) sites, labeled by j = 0, . . . M. 

Note, that the Cayley tree with q = 1 is a one dimensional lattice with N sites. 
The first non-trivial structure can be obtained for q = 2, which is described in the Fig. 
CD Here, in every level, each site is connected to two downward neighbours in the next 
level. There are M = (2 N — 1) sites which are labeled by j increasing from top to 
bottom and left to right. Further, each site of a Cayley tree can either be vacant or 
occupied by at most one particle. Correspondingly we define a variable Sj = 1,0 at 
site j, which represents a presence or absence of a particle respectively. These hardcore 
particles flow downwards (from level i to level i + with the following rules. A particle 
present in a level i, at the site say j, will move to its neighbour j' in the next lower 
level i + 1 only when Sj/ = 0. If more than one neighbour of j is vacant, one of them is 
chosen randomly for particle transfer. In the following, the bulk dynamics for q = 2 is 
described schematically. 
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Flow of particles is maintained in the system by the in- and out-fluxes of particles 
at boundaries. A particle can enter the system with rate a if the parent site j = 
is unoccupied (sj = 0). Any particle from the (q — 1) sites at the boundary level 
N — 1 can leave the system with rate 0. Note, that for q — 1, each site has only 
one downward neighbour and the dynamics is simply 10 — > 01. Such an exclusion 
process, namely Totally Asymmetric Simple Exclusion Process (TASEP)|8] on a one 
dimensional lattice, has been studied extensively for its application in vehicular traffic, 
cellular transport by motor-proteins, etc.. TASEP has been solved exactly [HI EE], where 
it shows novel properties like boundary driven phase transition, shock formation and 
propagation, condensation, and jamming. 

It is important to note that the bulk dynamics in the Cayley tree is particle 
conserving. In steady state, the average current that flow between neighbouring levels i 
and i + 1 is expected to be constant, say I q . Thus the average current on each link that 
connects a site at level % with its neighbour at level i + 1 is given by 

Ji^i+i = Ji = (3) 
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All configurations where a site in i th level is occupied and at least one of its (/neighbour 
in (i + l) th level is empty, contributes to the flow of particles. Thus, fraction of such 
configurations in steady state is qJi. Here, the factor q takes care of the fact that the 
average current I q is shared among q bonds. For example, when q = 2, we have 

2Ji = (l/\0) + + (0/\0), 

(4) 

where (. . .) represents the steady state averages. Similarly the total density at level i 



:*>■ (5) 



is the steady state average of the occupied sites belongs to level i. 

In the following we use the mean field theory (MFT) where both the fluctuations of 
particle density at the individual sites and variation of densities among sites in the same 
level, are neglected. At level i, thus, every site is assumed to have an average density 

Pi = —, 6 
Q 

where <pi is the total density at level i. Since a site of level i is occupied with probability 
Pi (and is vacant with probability 1 — p^), the average current through the bonds is given 
by, 

qJi = pi (1 - p? +1 ) , (7) 

where(l — p q i+ i) is the meanfield probability that at least one site of (i + l) th level is 
unoccupied. The factor q in front of Jj ensures that the outgoing current flows through 
q bonds. 

The conservation of particle density in the bulk of the system leads to a continuity 
equation 

-jjjT — J(i-1)-H ~ QJi-¥(i+l) — Ji-1 ~ £[Ji 

= ~A-i(1-pD-M1-P? + i)- (8) 
Q 

These equations must be supplemented by the following boundary conditions. First, at 
the root j = 0, where particle enters to the system, we have 

^ = a(l - pa) - qJ = a(l - p ) - p (l - pi). (9) 

Similarly at the last level i — N — 1, 

d PN-i T * 
— z — — J N-2 — PPn-i 
at 

= -pN-2^ - PaT_i) - 0PN-1- (10) 

From Eq. [H] it is clear that, in the steady state the bulk current is J, = Ji-i/q, 
which can be iterated to give 

J l = ~. (11) 
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Figure 2. Density profile for q — 2 obtained from Monte-Carlo simulations (symbols) 
with a — 0.7 and j3 = 0.2 are compared with Eq. [16] (dotted line) where po = 0.422 
was calculated using [T9l Clearly there are no finite size corrections. The densities of 
the last level pn-i also follow Eq. [20] (dashed line). In the inset we compare c versus 
/3 obtained from simulations (symbols) for a = 0.4, 0.7 with Eq. [2D] 



A comparison of this equation with Eq. ([3]) results, I q = qJo. Thus Eqs. (GO) and (J3]) 
provides an iterative equation for the density, 




First let us discuse the case q = 1 (i.e. TASEP) which is conceptually different 
from other cases q > 2 as the right hand side of the above map for q = 1 do not depend 
explicitly on i ( i.e. p i+ i = 1 — h/pi). In TASEP both, the density profile and the 
current show macroscopic changes as one vary the boundary rates a and 0, resulting in 
three different phases : (i) the high density phase(a > (3,(3 < |), (ii) the low density 
phase (a <!>,/? > a), and (iii) the maximal current phase (a > ~,(3 > |). TASEP 
is an exactly solvable model[_9j, however the meanfield analysis of Eq. [12] with q — 1 is 
known [TO] to produce correct phase diagram and the density profiles. 

One must treat the q > 2 case separately as the function in the right hand side of 
Eq. [121 changes in each iteration. The existence of the fixed points in such maps are not 
quite obvious. Let us take the large i limit and re- write the map in terms of (pi = q l pi, 
which is an 0{1) quantity. Thus, 




The fixed point of the map is then fa = 0*, given by 

fa = I q . (14) 

Thus in the bulk we have a constant level density fa = I q . This is similar to the case 
q = 1(TASEP), where the bulk density was found to be a constant. However, the 
density profile in TASEP show a boundary layer at both ends. Since the fixed point 
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solution <pi = I q is expected to hold only for large i, it is not clear if the total density 
4>i has a boundary layer for q > 2. The detailed mean field analysis (discussed below) 
show that the boundary layers are in fact absent for q > 2. 
To calculate I q we use the boundary condition (jSJ, 

J, = a(l-po). (15) 

Finally, Eqs. (]6]), flHJ) and ( IT5l) can be used to obtain the average level density, 



Pl = t = a{1 ' po) = h. (i6) 

qi gi qi 

Thus, the mean field densities p, and the current I q are expressed as a single parameter 
functions of p (which is same as cj) ). To calculate p we use the boundary condition 
OH]). It demands that, in the steady state p\ = i/l — Iq/ Po, whereas from Eq. [T6l we 
have pi = Ig/q. Hence, 



_ Q<1 ~ Po) = a(l - po) , 17 v 
Po q 

The above (q + l) th order equation in po can not have a close form solution for q > 3. 
However, numerical solution for any given a and q can be obtained with high accuracy. 
Analytical solutions of Eq. [T7| can be obtained for q — 2, 3. For example, when q = 2 
we have, 



2 V« 2 + 12a + 12 ,vr + # 

Po = - 1 cos{——)\, 18 

6 a 3 



where tcm^ = — — . (19) 



6V6 Va 4 + 18a 3 + 20a 2 + 24a + 8 
a 36 — 18a — a 2 

Thus, Eqs. (}T5]) and f|T6|) together with the solution of po from [19] provides the 
complete mean field solution of the asymmetric exclusion process on the Cayley tree. 
To check the validity of the MF theory, we simulate this exclusion process with a = 0.7 
and (3 = 0.2 on a Cayley tree (q = 2)with iV = 6, 8, 10, 12 and plot p, versus i in semi-log 
scale ( Fig. [2]). Clearly the densities p iy except pjv-i, show an exponential decay which 
agree quite well with the MF results f lT6|) . The plot depicts that the finite size corrections 
are absent here; i.e., the level density of, say at level 3 (i.e.p 3 ), is independent of the 
system size N = 6,8, 10, 12. Further, it appears that Pn-i varies with system size N as 
Pn-i = c/2 N ~ l , with c > po. 

To obtain the boundary density pn-i, we use the second boundary condition (ITU]) ; 
in steady state Jn-2 = Ppn-i- Since Jn-2 = Iqlq N ~ x [from Eq. |3], we have 

^-i = qFn Wlth c = J = p • ( 2 °) 

In the inset of Fig. H] c versus /3 obtained from numerical simulations for two different 
values of a = 0.4, 0.7 are compared with corresponding mean field values given by Eq. 
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Figure 3. Plot of J q verses a with different q = 1,2,3,4.. For q = 1 (TASEP) the 
current J\ depends on (3 which is taken to be j3 = 0.2 Inset : For q = 2, J versus a 
obtained from numerical simulations with j3 = 0.2, 0.6 (symbols) is compared with Eq. 

ma 



For a generic Cayley tree with q > 2, the current I q = a(l — po) where p is given 
by Eq. ( JT9l) . Evidently, I q is independent of the exit rate /3 and increases monotonically 
with a. To verify this we calculate I q for a Cayley tree with q = 2 using Monte-Carlo 
simulations for two different values of (3 = 0.2, 0.6 and plot them against a (inset of 
Fig. [3]). The mean field current (fl5|) for q = 2 , draws as a solid line there, shows an 
excellent agreement. Current I q for q > 2 are shown in the main figure [results from the 
simulations data (not shown here) matches very well with Eq.l TToT) ]. 

Figure [3] compares I q for different q = 1,2, 3, 4. It is only Ii, the current in TASEP, 
that depends on (3 (chosen here as (3 = 0.2). As expected, I q is a strictly non- decreasing 
function of the entry rate a. Thus, the maximum achievable current on a Cayley tree 
is I q (a = 1). For example, when q — 2, the maximum current 

ir x = 3 [1 + Wcos( - + -tan~\^-) )] = 0.485 (21) 

is almost twice as large as that of TASEP (q = From Fig. [31 it appears that 

lim^oo I q (a — 1) = |. This can be understood from the fact that when q — >■ oo, the 
rate of out flow from the root j = is unity, as one of the the infinitely many neighbours 
of the root are expected to be vacant with probability 1. Thus the particle density p 
at the root is expected to be 1/2 for a = 1, resulting in 1^ = 1/2. Another interesting 
fact about I q , is the following in-equality, 

Ji(a,/3)</ 2 (a)...</ 0O (a). (22) 

which holds for any arbitrary value of (a, (3). 

In conclusion, we have studied the asymmetric exclusion process on a regular Cayley 
tree with arbitrary coordination number z = q + 1, where particles enter the system 
only at the parent site with rate a, and exit from one of the sites at the last level with 
rate (3. In the bulk they move from occupied sites to one of their downward unoccupied 
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neighbours chosen randomly. TASEP is a special case of this model for q — 1 which 
exhibit boundary driven phase transitions. In this case the steady state current I q is 
different in different phases and depend both on a and (5. Surprisingly, for q > 2, we find 
that there is only one phase and current I q is independent of (3. Further, I q increases 
monotonically with a and reaches a maximum value twice as large as that of TASEP 
even for small coordination numbers (say, q — 2). Again, the density profiles do not 
show any boundary layers or finite size corrections. It is like a free flow of particles from 
one end to the other; only the last level is effected by the exit rate. The model could 
find application in nutrient transport in blood-circulation systems, air circulation in 
lung or in antibody transport in immune systems or water transport in trees. Possibly, 
branching is essential in all these systems (say trees) as the flow (of water at different 
levels) need to be maintained independent of the exit rates (weather conditions) at the 
last level (leaves). 
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